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Abstract 

'"17' The rotation-less Newton-Hooke - type symmetry found recently in the Hill problem and instru- 

' mental for explaining the center-of-mass decomposition is generalized to an arbitrary anisotropic 

(— I oscillator in the plane. Conversely, the latter system is shown, by the orbit method, to be the most 



general one with such a symmetry. Full Newton-Hooke symmetry is recovered in the isotropic case. 
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1. INTRODUCTION 



Renewed interest in Kohn's theorem on decomposing a system of charged particles in 
a magnetic field into center-of-mass and relative coordinates stems from relating it to the 
Newton- Hooke symmetry of the Landau problem P-Ej. Recently, another example was 
found, however ^ |7j : Hill's equations studied in Celestial Mechanics also admit a center- 
of-mass decomposition, but no full Newton-Hooke symmetry. Rotations are broken, but 
translations and (generalized) boosts are still symmetries, hinting at that it is the sub- 
group spanned by the latter which is important for the purpose; additional symmetries like 
rotations are irrelevant. 

In this paper, when referring to Newton-Hooke type symmetry it is Newton-Hooke with 
or without rotations that we have in mind. 

The possibility of decomposing an isolated system into center-of-mass and relative co- 
ordinates has been linked to Galilean symmetry : Souriau [8] argued, in fact, that this 
property depends on the Galilei group having an invariant Abelian subgroup — namely the 
one spanned by translations and (Galilean) boosts. 

Remarkably, the cohomological structure which determines the existence of a central 
extension originates precisely in this subgroup |8]. Remember that, in dimension d > 3, 
both the Newton-Hooke and the Galilei groups have a one-parameter central extension, but 
in the plane they both admit an "exotic", two-parameter central extension [SHI2]- 

In this paper, we focus our attention at the Newton-Hooke case, the Galilean one being 
rather standard. 

The ordinary Landau problem admits the one-parameter centrally extended Newton- 
Hooke group as symmetry [Ij, but the "exotic" [non- commutative] version has indeed the 
two-parameter version |;4j. In the Hill case, rotation-less "Newton-Hooke type" symmetry, 
with one (or, in the "exotic" case, with two) central extensions could be established [HI Ej. 

Our main result proved in Section [5} is : 

Theorem 1 : The most general planar system with Newton-Hooke - type symmetry is a 
(possibly non- commutative) anisotropic oscillator in a uniform magnetic background. The 
symmetry extends to full Newton-Hooke symmetry in the isotropic case. 

The proof will be accomplished by applying the orbit method, which provides us indeed 
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with all systems upon which the symmetry group acts transitively [H [T3HT5] . 

From the technical point of view, we will find it convenient to use chiral decomposition 
P m [121 [ISl [13 US\ ■ The motion in the (ordinary) Hall effect can in fact be decomposed into 
two uncoupled chiral oscillators with opposite chirality [i6\. Conversely, combining two Id 
chiral oscillators may yield the non-commutative Landau problem [T7H20] : then the chiral 
method allows for an elegant derivation of the (Newton-Hooke) symmetry. 

Recently, the method was extended to the Hill problem |6l [7] which is in fact a "maximally 
anisotropic oscillator" ; here we further extend it to arbitrary anisotropy. 

Our paper is organized as follows. 

In section [2] chiral oscillators are reviewed. 

Then we outline the Newton-Hooke symmetry in the Landau- , and the rotation-less 
Newton-Hooke type symmetry in the Hill cases, respectively. 

In Section |4| we generalize to an arbitrary, possibly anisotropic, oscillator. 

In Section [5] we proceed conversely: applying the orbit method we describe all systems 
with Newton-Hooke-type symmetry acting transitively. 

We also study the arising of further symmetries and explain the difference between sym- 
metry with or without rotations. Our results allow us to deduce 



Theorem 2 : The system is either a truly anisotropic oscillator with Newton-Hooke-type sym- 
metry only and no rotations, or it is isotropic with full Newton-Hooke symmetry, including 
rotations. 

Moreover, in the first case, it can be brought into a "Hill-type form" , and in the second 



one it can be transformed into a free particle, cf. Sections 5 D and |6] 
An outlook is presented in the Conclusion, section [7} 



2. CHIRAL OSCILLATORS 



Chiral oscillators arise owing to the ambiguity of the phase-space description of a har- 
monic oscillator. In detail, let us consider a one-dimensional harmonic oscillator of unit mass 
m = 1 and frequency u. Viewing the position and velocity, x and i;, simply as coordinates 
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FIG. 1: The phase-space trajectory of a chiral oscillator turns clockwise or anti-clockwise, depend- 
ing on the sign of the frequency. Both trajectories project, however, onto the same motion in 
configuration space. 



on the phase space, we write 

yi = x, y2 = i, 
and consider the two first-order phase space Lagrangians 

The associated (Euler-Lagrange) equations read 



(2.1) 



(2.2) 



(2.3) 



Our clue is that, for both signs in eqn. (2.2), eliminating either yi or y2 yields, for the 



remaining variable, the same equation, namely that of a Id harmonic oscillator. 



jji + u yi = 0, 



i = 1,2. 



(2.4) 



The solutions of (2.3) are simple rotations in phase space - but in opposite directions, 



depending on the sign [28j. (This is indeed the very meaning of the word "chiral"). Then 
we note that both types of motions project into configuration space according to the same 
motion x{t), as illustrated on Fig. [l] 
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We note that the same conclusion can be reached using a Hamiltonian framework. The 



eqns. (2.3) are indeed those of the symplectic structure and Hamiltonian 



n± = ±\eijdyi A dy 



-1 01 



H = luy' 



(2.5) 



namely iji = [yi,H^^, where the Poisson brackets { ■ , -^^ are those associated with the 
chosen symplectic structure Q±. The coordinates yi are non-commuting, 



{yi,y2}_ 



(2.6) 



as it is natural for position and momentum on the phase space. We mention for com- 



pleteness that the Lagrangians (2.2) are the Cartan forms of the Souriau forms [HI 1^ . 



L^dt = A 



dX+ = n+-dH A dt. 



(2.7) 



3. TWO SYSTEMS WITH NEWTON-HOOKE [TYPE] SYMMETRIES 
A. Landau problem 

The classical example of a system with one-parameter-centrally-extended Newton-Hooke 
symmetry is provided by the "ordinary" [meaning commutative] Landau problem [1]. Gen- 
eralizing the latter, we consider "exotic" particles endowed with masses, charges and 
non-commutative parameters rria, and 6a, respectively, moving in a planar electromag- 
netic field B,E Following [HI [T2], we describe our system by 

(3.1) 

Pi = CaBe'^xi + eaE\ 

where 

m* = A^ma with = 1 - CaOaB (3.2) 

is the effective mass of the particle labeled by a = 1, . . . , A^. Note, in the first relations, 
also the "anomalous velocity terms" perpendicular to the electric field, E. The variables Pa 
here could be called "momenta" - but to avoid confusion with the conserved quantities, we 
simply consider them as coordinates on the phase space. 

Although our theory works for any B and E, we assume, for simplicity, that the magnetic 
field is constant, B = const, and that the electric field is that of an isotropic harmonic 
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trap, —kx, augmented with an interparticle force coming from some two-body potential, 
For 6a = 0, the ordinary Landau problem is recovered. 

Summing over all particles, we find that when ea/rria and CaOa are both constants i.e. 
when the generalized Kohn conditions [3], 

t^^ = ^ = ^^ = n, eaea = ee, e = ^",7."^" (s-s) 

ma M ' ' M2 ^ ^ 

hold (where M = J2a ^a. and e = J2a ^a. are the total mass and charge), then the center-of- 
mass, X = '^^niaXa/M, splits off, 

(3.4) 

pi = eBe'^X^ + eE\ 

where 

M* = AM, A = 1 - eQB, P = ^pa. (3.5) 

a 

The center-of-mass behaves hence as a single "exotic" particle carrying the total mass, charge 
and non-commutative parameter, M, e and B, respectively. 



We note that eqns. (3.4) are in fact Hamilton's equations for the Poisson structure [19], 
{X\X^} = -e'\ {X\P^} = -, {P\P^} = —e^K (3.7) 



When the [generalized] Kohn condition (3.3) is satisfied, then, for identical initial veloc- 
ities, all individual particles move in the same way, — and this motion is shared by their 
center of mass, cf. Fig. |2} 

The best way to understand the intuitive content of the Kohn condition is, however, to 
consider what happens when the Kohn condition is not satisfied. Consider, for example, 
two particles in a pure magnetic field such that 

= — = 2ki = 2 — . (3.8) 

m2 nil 

Then, assuming identical initial velocities, each of them performs a rotational motion but 
with different radii, 

R={m/e)^ => R2 = IRi, (3.9) 
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FIG. 2: If the Kohn conditions (3.3) are satisfied, all particles turn along circles of equal radii with 



common angular velocity. Their motion is shared by their center of mass (in dashed). 
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FIG. 3: // the Kohn conditions (3.3) are not satisfied, K2 = e2/m2 = 2ei/mi = 2ki, for example, 
the individual radii and the angular velocities are different. The motion is not more collective, and 
the center of mass describes a complicated (dashed) curve. 



and with different frequencies: 



V 

R 



U2 = 2ui 



(3.10) 



[so that uiRi = V = U2R2]- Their center-of-mass would then clearly not move on a circle, 
rather on some complicated curve, cf. Fig. [3j The 3-body situation is illustrated on Fig. |4} 



Symmetries . 

Let us restrict ourselves henceforth to the purely magnetic case, E = and to the 
center-of-mass motion. The coordinate P is not conserved; one readily shows, however, 
that the "magnetic momentum" (which can also be derived by Noether's theorem as the 
conserved quantity associated with the translational symmetry [20]) [21] and "magnetic 
center-of-mass" , 

(3.11) 

K = MA^R{uj*t)X, 
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FIG. 4: The behavior of a 3-hody system, (a) If the Kohn conditions (3.3) are satisfied, all particles 



move collectively, along with their center- of -mass, (h) If (3.3) is not satisfied, ki = 1, K2 = 2, = 
3, for example, the motion is not more collective, and the center of mass (in dashed) describes a 
complicated curve. 

respectively, where u* = eB/M* = u/A, are both conserved [30], and span indeed two 
uncoupled Heisenberg algebras with central charges —Mu and AMu, 

{U\ W} = -Mue'^, {K\ K^} = AMuje'\ K^} = 0. (3.12) 

Time translations and rotations are plainly symmetries also, and the associated conserved 



quantities, namely the Hamiltonian H in (3.6) [3T], augmented with the total angular mo- 
mentum [m [in] , 

J = XxP + ^x' + ^P\ (3.13) 

have commutation relations 

{E, n'} =0, {if, K^} = e^^K\ (3.14) 

{J,X^} = £*^X^ {j,P')=e'^P\ {J,i7}=0. (3.15) 

In conclusion, the exotic Landau problem [with our without an isotropic harmonic trap- 
ping force] admits an "exotic" i.e. two-parameter centrally extended Newton-Hooke sym- 
metry in [20]. In the commutative case G = 0, the central charges are correlated, =fMci;, 
and the symmetry reduces to the one-parameter extension studied in [1]. 



We record for further use that the total angular momentum, J in (3.13), can also be 



presented in a number of different ways. Firstly, we note that the new variables [TT] 

Qi = Xi + ^{l~ Vl - deB^ Eijpj , (3.16) 

(3.17) 



1 + VI - OeB 

Pi = Pi 



are canonical, and in their terms the total angular momentum is simply 



J = Qx P. 



(3.18) 



Here we just mention that using chiral coordinates (sect. 4C), the angular momentum can 



further be decomposed, see (4.22). 



From now on the generalized Kohn conditions (3.3) will always be assumed, allowing 



us to consider the center-of-mass alone. Coordinates will again be denoted by lower-case 
letters, as for a one-particle theory. 



B. Hill problem 

The center-of-mass of celestial bodies, moving approximately along circular orbits in 
the gravitational field of some heavy central object, is described, in the linear approximation, 
by the planar equations [B], 

X — 2uy — Su'^x = 0, 

(3.19) 

y + 2ux = 0, 

where u is the angular velocity = GM/R^ for a circular Keplerian trajectory of radius 
R. Note here the anisotropic oscillator term, which is the remnant of the centrifugal and 
Newtonian expressions under the used approximation. 

The system is conveniently analyzed in terms of the chiral coordinates p] 

x = X+ + X^, = ^u;Xl, p'^ = -2uXl - '^ujX^-- (3.20) 

Our investigations can in fact be generalized to exotic particles with 6* 7^ 0, see [7]. Skipping 
details, we mention that in both [exotic or not] cases we find that ordinary translations and 
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n 



certain "time dependent translations" (also called "generalized boosts"), 

^ Xl{t) cos{uj/A)t - ^Xl{t) sin(w/A)t ^ 
y 2TXl{t) sm{uj/A)t + Xl{t) cos{uj/A)t j 



(3.21) 



K 



are conserved, where 



A = l-2mu;^, T = l- 3m9uj/2. 



(3.22) 



Their commutation relations are, once again, those of two exotic Heisenberg algebras with 
central charges —{2/muj) and (r/A)(2ma;), respectively. 



moo 



{K\K'} 



r 2 

A muj 



{n„ir,} = o 



(3.23) 



In the commutative case F = A = 1, and the one-parameter centrally extended symmetry 
found in [6] is recovered. 
The Hamiltonian, 



H = H.+H_ 



2 ■ 4r2' 

is also conserved. Its commutation relations with translations and boosts read 



(3.24) 



{H,U'} = 0, 
{H,K^} - 



2f 



(3.25) 



K'\ {H,K^} = 2ru*K^ 



As rotational symmetry is plainly broken, the total symmetry of the Hill problem is 
Newton-Hooke without rotations. 

4. ANISOTROPIC HARMONIC OSCILLATOR 



We note that the Hill problem is in fact a maximally anisotropic "one sided" oscillator. 
The case of a general anisotropic oscillator is worth studying in some detail therefore. 
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A. Chiral coordinates 



Consistently with the general theory sketched in Section 3 A, an "exotic" [i.e., non- 
commutative] charged harmonic oscillator in the plane in a homogenous magnetic field B is 
described by the symplectic form and Hamiltonian, 

Q = dp' A dx' + -e'Hp' A dp> + —e'^dx' A dx^ , (4.1) 



„2 

H = ^ + V, V = \kixl + \k2xl, (4.2) 



respectively, with the parameters having the same physical interpretation as before. The 
spring constants ki and ^2 may or may not be identical. 

The idea of Alvarez et al. [17] has been to combine chiral oscillators. Multiplying both 
the symplectic form and the Hamiltonian (or, alternatively, the Lagrangian) by the same 
overall constant /x, 

fi— i-yui^, H-^fiH i.e. L^fiL, 

would not change the equations of motion. But what happens, if we multiply them with 
different coefficients before adding them ? Conversely, can we decompose a given system 
into two chiral parts ? To answer these questions we introduce, following [6l [T71 [20], new 
coordinates on the phase space. 



X 



X; + Xi, p^ = a+Xl + a^X^, p'^ = -f3+Xl - f3_X^_, (4.3) 



where the coefficients a± and P± will be determined from the requirement that both the 
symplectic form and the Hamiltonian should split into two uncoupled one- dimensional sub- 



systems we shall call chiral components. Inserting (4.3) into (4.1) shows that the symplectic 
form splits as 1] = + , whenever 

a- + /3+- ea-(3+ = eB, a+ + /3_ - 6'a+/3_ = eB. (4.4) 



Similarly, inserting (4.3) into (4.2) yields that the Hamiltonian splits as if = H^ + H^ when 



+ mk2 = 0, + fnki = 0. (4.5) 
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Then a tedious calculation allows choosing 

1 

'2{eB + 9mki) 



a+ = , A I -e^B^ + m{k2-ki) + e^m''hh (4.6) 



+ \jAmki{eB + 9mk2)'^ + (e^B^ - m{ki - ^2) - 9'^m'^kik2Y^ , 



and 



+'\J 4mki{eB + 9mk2Y + (e^S^ - m{ki - k2) - 9'^m?kik2)^ , 



4mki{eB + 9mk2Y + (e^S^ - m{ki - k2) - 9'^m?kik2)' 
provides us with decomposed symplectic form and the Hamiltonian, 

Vt = + (4.10) 

{-a+ - /3+ + 9a+l3+ + eB) dXl A dXl + (-a_ - /3_ + + eB) dX]_ A dX! , 



— V 



and 

= H+ + H_^ — X (4.11) 
2m 

(/32 + mfci) XlXl + (a^ + ^^Js) + (/^^ + ^A;i) X^_X^_ + (a^ + ^A;2) X^_X^_ 

respectively. 
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For ^ = the commutative cases [Gj [T7] are recovered; 

For ki = —3mu'^, k2 = and B = 2u, we get 

moj 3 
a+ = 0, a_ = — , (3+ = -mu, f3_ = 2ma;, 

(4.12) 

mu A muj 



2 ' ^' 12' 
and the results found before in the Hill Problem [HI E] are obtained; 



When ki = k2 our oscillator is isotropic. Then a± = /3±, and (4.10 4.11) reduce to the 



chiral decomposition for the [exotic] Landau problem with harmonic force, studied in 



For ki = k2 = the oscillator is switched off, and the system reduces to the purely- 
magnetic non- commutative Landau problem fTli [T^ fIU\. 



B. Motions 



Let us assume that none of the coefficients fi± vanishes. Then it follows from (4.11) that 
our chiral coordinates satisfy the Poisson bracket relations 

{Xl,Xi] = -—e'^, {Xl,Xi]=0 {Xt,Xi] = -—e'^ . (4.13) 

The equations of motion read therefore 

m/i-tX4 = — (a|+mfc2)X^, 

m/x±X| = {l3\ + mki)X\. 

Both chiral components X± are governed, hence, by uncoupled equations which are rem- 
iniscent of those of Id harmonic oscillators, to which they reduce, however, only in the 
isotropic case, ki = k2. 



(4.14) 



Assuming a\ + mk2 ^ [32], eqns. (4.14) are solved by 
X]. = A± cos co±t + B± sin u±t, 

X| = F± (^A± sin uj±t — B± cos CO ±tj, 



a% + mk2 ' 
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FIG. 5: "Epicyclic" motion of an oscillator in a constant magnetic field, (a) is the isotropic case 
ki = ^2; <ind (b) is an anisotropic case with ki = 4k2. The black line is the physical trajectory, 
X = + X_, and the dotted red line is the motion of the guiding center, X^{t). 

where the frequencies read 



m/i-t 

Both X-t-trajectories are elhpses, as illustrated in Fig. |5} Note that the frequencies, u+ and 
are in general different even in the isotropic case, and the curves do not close therefore. 



C. Symmetries 



Eqns. (4.15) allow us to infer that 



1 



A± = X]_cosoj±t+—X_^smco±t, 

1 1 2 

B± = s'mu±t — -^X^ cos u±t 
F± 



(4.17) 



are conserved. A direct calculation yields, furthermore, for both labels ±, the uncoupled 
Heisenberg algebra relations 

1 



(■)- 



0. 



(4.18) 



Adding the Hamiltonian (4.11), the doubly-centrally-extended rotation-less Newton-Hooke 
algebra is obtained. 



15 



Both sets of chiral coordinates X± describe 2d symplectic vectorspaces. The symplectic 
forms VL± are plainly symmetric under phase-space chiral rotations, X± — )■ R{X±). None of 
the Hamiltonians H± is symmetric in general, though. The natural diagonal action. 



x = X+ + X_^R{X+) + R{X_) = R{x), (4.19) 

is not a symmetry therefore : rotations are broken by the anisotropy. 
In the isotropic case, 

ki = k2, (4.20) 

however, we have a± = P± and the coefficients of the quadratic terms both in and if_ 
are hence identical, so that the chiral rotations X± — )■ R{X±) do act as symmetries for 



the components : rotational symmetry is restored. The square-root factors in (4.15) become 
unity, F± = 1, and the trajectories become circles. The frequencies, 

ai + mk 

oo± = , 4.21 



are not identical, though, since a+ 7^ a_ and /i+ 7^ in general, cf. (4.6) - (4.7) and 



(4.10) 



It is worth recording that, in terms of chiral coordinates, the total angular momentum. 



(3.13), is also decomposed, as 



J=J+ + J_, J+ = _(^X+j , J_ = -A— (X_j , (4.22) 

where A = 1 — cBd, as before. Its conservation, J = 0, can also be checked directly, using 
the equations of motion. 

We just mention that the singular case = or yU_ = 0, leading to Hall-type motion, 
can be dealt with as in the previous occasions, [U U\ [HI ttH 120] . 



5. SYSTEMS WITH PRESCRIBED NH-TYPE SYMMETRY 

We wish to show now that the results of the previous Sections fit, in fact, into a general 
framework. To this end, we assume that the dynamics under consideration is invariant under 
the transitive action of some Lie group G, and then classify all such symplectic manifolds 
upon which G acts by symplectically. The proper tool for doing this is provided by the orbit 
method ^^EM- 

16 



Our choice for the group G is dictated by the following considerations. As far as possible, 
we would like to allow for a generalization which includes both the Galilei and the Newton- 
Hooke groups, and also the "rotation-free part" of the latter, considered in Sect. [3j The 
main characteristic features are therefore the following: 

(i) there exists generators (namely of boosts and momenta) which, via the orbit method, 
yield the basic canonical variables; 

(ii) The Hamiltonian equations of motion are linear in the latter variables; the Hamilto- 
nian belongs therefore to the Lie algebra itself, and acts linearly on the remaining variables. 

We want our generalization to be a minimal one in that no further symmetry generators 
beyond the above ones should be included. Such generators will appear later however for 
specific values of the structure constants. 

Guided by these considerations, we start with the following Lie algebra commutation 
relations. 



= lu.jM, i,j = l,...,2N 

[M,e.] = 0, 



(5.1) 



[M,H] = 0, 

where u = (wjj) is a non-singular antisymmetric matrix. The only non-trivial Jacobi identity, 

(cyclic) = 0, (5.2) 
yields the constraint Ai^Ukj — AjkUj^i = 0, i.e., that B = Au is a symmetric matrix, = B. 



The algebra (5.1) admits the Casimir operator of the form 

C = MH- (5.3) 

where without loss of generality we can assume that X = (Xij) is symmetric. C commutes 
with all generators, provided A = —uX. 
Collecting our results, our algebra reads 

[M,(-)] = 0, 

(5.4) 

= —iuJikXkj^j, 
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and is uniquely defined by choosing the non-singular antisymmetric matrix oj and the sym- 
metric matrix X. 



The next step is to classify the inequivalent algebras (5.4). Under the invertible transfor- 
mation 

= det (A,) 7^ 0, (5.5) 

The matrices u and X transform according to 

uj' = DuD'^, X' = {D-^fXD-^ . (5.6) 



Using the latter we can find the "canonical" form in any class of equivalent algebras (5.4). 

In what follows we shall restrict ourselves to the case 2N = 4, the generalization to 
arbitrary being straightforward. 

To complete our classification scheme some further assumptions on the matrix X have 
to be made. The existence of the Casimir operator C implies that, on each orbit, the 
Hamiltonian is a quadratic function of the basic canonical variables, to which a trivial term, 
representing the internal energy, has been added (see Appendix B). Whether the energy is 
positive definite or not depends, therefore, on the choice of X. 

The following cases will be considered separately. 

A. X Positive definite 

Consider first the case of a positive definite matrix X. By an appropriate choice of D in 



eqns. (5.6), X = I can be achieved. In fact, X, being symmetric, can be diagonalized by a 
suitable orthogonal transformation. Then an additional diagonal transformation reduces X 
to the unit matrix. 

Assuming X = J, we still have some residual transformations left at our disposition. 



Namely, as it is seen from eqns. (5.6), D can be taken to be an arbitrary orthogonal matrix, 
without spoiling the condition X = I. The question is now to classify all antisymmetric 
4x4 matrices u up to an orthogonal transformation. This problem is solved in Appendix 
A (which is actually the Euclidean version of the classification problem for electromagnetic 
field configurations under the action of the Lorentz group). As shown in Appendix A, u can 
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be put into the form 



UJ 





-VLi 
V -^2 



Vl2 







fil2 > 0. 



Defining 



one finds the following non-trivial commutators : 

[Bi, Pk] = iSikM, 
[H, Bi] = -zP„ 
[H,P,]=tQ]B„ 



together with 



C = MH - \ [Pi + P| + ^\B\ + ^\Bl) . 



(5.7) 



(5.1 



(5.9) 



(5.10) 



Orbits 

We can now apply the orbit method (Appendix B). Consider the coadjoint orbit 
parametrized by m > 0, the coordinate in dual space corresponding to the Casimir op- 
erator M and interpreted as the mass, and by em, corresponding to the Casimir operator C 
and interpreted as the internal energy. Let Pi, bi, h, i = 1,2 he the relevant coordinates in 



the space dual to the Lie algebra (5.9 ). As shown in Appendix B, the points of the coadjoint 



orbit are parametrized by pi and bi . Defining 



qi = bi/m. 



(5.11) 



we find 



p{_ 
2m 



P2_ 

2m 



(5,12) 



Hence, we arrive at an in general anisotropic oscillator, as the most general system with 
positive definite energy, admitting the symmetry defined by the rotation-less Newton-Hooke 



commutation relations (5.1). 
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B. X semi-positive 



Let us consider the case when the matrix X is semidefinite. We restrict ourselves to X 
having a single zero eigenvalue (as in the Hill case). Then one can select the matrix D in 



(5.6) in such a way that X acquires the form 

/ 1 





X 







v 





1 
1 / 



(5.13) 



One can show again (see Appendix A) that the residual freedom in the choice of the basis of 



our algebra allows us to put u into the form (5.7). Using again eqns (5.8), we find therefore 



[H, Pi] 
C 



(5.14) 



MH - i (P2 + p2 ^ ^2^2 j 

The orbit method yields, in this case, the dynamics describing a harmonic oscillator in one 
direction, and free motion in the second one — as in the Hill problem [6]. 
The case of multiple null eigenvalues of X can be dealt with similarly. 



X indefinite 



Let us drop, finally, the assumption of positive (semi)definiteness of X. We consider in 
more detail the cases of two positive - one negative - one null eigenvalues. By an appropriate 
choice of D one can achieve 



X 




where G = diag(— 1, 1, 1). 



(5.15) 



According to the results in Appendix A, the symplectic form u can acquire three canonical 



forms, namely those presented in eqns (A. 5) - (A. 16) - (A. 17). Then the orbit method gives 
the following dynamical systems : 
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(ii) 



{(li,Pk} = Sik, (5.16) 



{qi,qj} = aeij, {pi,pj} = reij, (5.18) 
h= ^+(^-'i)+e; (5.19) 



111 



{qi,qj} = aeij, {^2,^2} = 1, {Pi,Pj} = retj, (5.20) 

2m 2my ^ ^ 

The parameter a here is a clear indication of n on- commuting nature of the coordinates qi 
and g2 [S3]- 

The case of non (semi) definite Hamiltonian is the most involved one. Unlike in the pre- 
vious cases, after the "canonical" Hamiltonian is fixed, there still remain three inequivalent 
forms of the basic Poisson brackets. 

The reason for that is clearly seen from the derivation given in Appendix A. The 3x3 
submatrix Ug of the matrix u transforms, under the transformations leaving the form of 
the Hamiltonian invariant, as an 0(2, 1) antisymmetric tensor. Its canonical form depends 
therefore on the value of the "electromagnetic" invariant 

2 

^(woi)' - (W12)'. (5.22) 

i=l 

Depending on its value, the basic Poisson brackets can take different, inequivalent forms 



(assuming the form of Hamiltonian is fixed). The labeling of variables in equations (5.16) 



(5.21) is dictated by our preference for the form of the Hamiltonian, rather then that of the 
Poisson brackets. It must be stressed, however, that the final choice of appropriate variables 
should be dictated by additional assumptions, not resulting from symmetry considerations 
only. 
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As an example, let us consider the planar Hill equations, as presented in Refs. [6l [7] . The 
Hamiltonian reads 

^=^(P?+P2)-— 92, (5.23) 
and yields Hill's equations for the following Poisson brackets, 

[Qi^Pj] = ^ij: [Pi^Pj] = '^mueij , (5.24) 

where the commutative case, a = 0, has been chosen for simplicity. The other parameter is 
= Sm^w^, and B = 2muj is the effective "magnetic" field. Let us put 

6 = Agi, C.2 = VSmuj q2, ^3 = ^, U = -^ (5.25) 



with A 7^ arbitrary. Then H acquires the standard form 

ii = \{e,^ii-e^. (5.26) 

and the relevant Poisson brackets read 

{6,^4} = V3c^, {i2.U) = 2cc;. (5.27) 
Therefore one finds, with the notations of Appendix A, 

Uoi = 0, Wo2 = VSUJ, Ui2 = 2u, — bj\2 = —uj"^ < 0. (5.28) 



According to the classification given in Appendix A, we are dealing with the case (A. 14), 



and the "canonical" form of the Poisson brackets is given by eqn. (A. 16), in full agreement 
with the results of Refs. [6],[7]. 

D. Additional symmetries 

We now study the question of additional symmetries. Consider the case of a positive 
definite Hamiltonian. As it has been shown in the previous Section, the initial algebra can 
be put into the form 





(5.29) 




(5.30) 


[M, .] = 0, 


(5.31) 
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where u is given by eqn. (5.7). We add a new generator J which is assumed to obey 



[J,M] = 0, [J,if] = 0, [J,^i\=ij^k^k. 
where j = (jik) is an appropriate matrix. The two additional Jacobi identities 



(5.32) 



[J, [H, e.]] + (cycl) = [e„ [J, 0]] + (cycl) = (5.33) 
yield ju + uj'^ = 0, ju — uj = 0. Hence j = —j'^, and the general solution reads 
(i) ^1 ^ ^2, 



( 



] = a 







-1 



1 \ 




\ 



+ /3 







\ 

1 





V -1 



(5.34) 



(ii) Vli = 
( 



] = a 



1 \ 







-1 
\ 







-/3 





V -1 



\ 



1 







+7 



1 \ 



1 



-1 

V -1 







/ 1 



+5 



J 



\ 



-1 







1 



V 



-1 / 
(5.35) 

Before explaining the meaning of the particular solutions, let us note that, once the 



equations (5.33) are obeyed, there exists a second Casimir operator, namely 



(5.36) 



where Y = —u ^j. If the coadjoint orbit is parametrized by the value ma of the Casimir C, 



eqn. (5.36) |[3l] yields the expression for J which, as in the case of the Hamiltonian, consists 



of the sum of the quadratic term plus the "internal" contribution to J, 



mJ = lYijZiZj + ma, 



(5.37) 



where the z^s are the basic variables parametrizing the points on the orbit (cf. (B.l) in 



Appendix B). Eqn. (5.37) allows us interpret a as the internal angular momentum, analogous 
to the internal energy, e, in our previous formulae, cf. [8]. Using the general solution for 
j, eqns. (5.34) and ( 5.35[ ), one finds the following generators as functions on the coadjoint 
orbit (up to an internal part) : 
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(i) for VLi^VL2, 



i^i-^ + ^gi, H,-— + —q,, (5.38) 

(ii) for Qi = Q2 = ^, 

pj mfi2 pI ^^2 

-"1 = ;S \ T^Qli ^2 = 7) \ 7^Q2^ 

2m 2 2m 2 gg-j 

J =qiP2-q2Pi, Z = pip2 + m^Q?qiq2. 

The meaning of the above expressions is clear. First of all, for an anisotropic oscillator 
we have two integrals, corresponding to the partial energies; our system is integrable. 

For equal frequencies, the dynamics is superintegrable : there are three functionally 
independent integrals. One can choose the angular momentum as the third one. The 



four integrals in eqn. (5.39) are linearly independent but they are functionally dependent. 



Note also that our integrals (5.39), being quadratic in canonical variables, form the u(2) 
Lie algebra — the well-known dynamical algebra of a two-dimensional isotropic harmonic 
oscillator. In fact, if one defines 

V2^^{H2-H^), V,^^ip,p2 + m'n'q,q2), (5.40) 

the resulting Poisson brackets algebra reads 

{Vi,Vj} = eiji^Vk, (5.41) 

i.e., span the su(2) algebra. The fourth generator, namely the Hamiltonian, 

Vo = H, (5.42) 

can also be added [I?]. Vq commutes with all other V^s, completing the su(2) algebra into 
the unitary algebra u(2). 

Let us remark that even for Qi 7^ Q2 there exists an additional integral, provided the 
ratio of the frequencies is rational, 

r = fii/fis = m/n. (5.43) 

It is, however, no longer quadratic in the canonical variables, yielding a W^-algebra, instead 
of a Lie algebra [22j. In fact, an additional integral of the motion which yields our system 
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superintegrable can be constructed as follows. One defines the classical counterparts of the 
creation/annihilation operators by 

0-1 = qi 7^, ai = qi^ — . (5.44) 

It is then easy to check that 

C"'"^ = (ai)"(a2)™ (5.45) 
is an integral of the motion. In the isotropic case n = m = 1, for example, 



is a combination of those conserved quantities in the second line of (5.39), namely of the 
angular momentum and the "mixed" quantity denoted by Z. 

The integral C**^'™ is functionally independent of the partial energies. Hi 2- Moreover, 
there are no further independent (and explicitly time-independent) integrals; therefore, the 
Poisson bracket between H12 and C*"'™ are functionally expressible in terms of them, and 
form a finite ly-algebra 



Let us conclude this section with some remarks. We have shown, at least in the case of 
(semi) definite hamiltonian, that there exists a unique "canonical" form of the underlying 
dynamics. However, the choice of this canonical form is dictated by mathematical simplicity 
rather than by physical requirements which are, in fact, additional assumptions. It seems 
reasonable to assume, generally, that the physical variables are those which convert the 
system into (non-commutative) anisotropic oscillator in a uniform magnetic background. 
This can be always done because our canonical form may be converted back into any other 
hamiltonian form obeying the symmetry assumptions. Therefore, we end up with Theorem 
1, as stated in the Introduction. 



6. THE BARGMANN POINT OF VIEW 



The NH symmetry of an isotropic oscillator can conveniently be derived by "importing" 
the Galilei symmetry of a free particle using Niederer's transformation, which maps every 
half period of the oscillator onto a free particle [3l Elj. One way of seeing this is to work 
within Duval's "Bargmann" framework, where classical non-relativistic motions are null 
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geodesies of a suitable relativistic spacetime [251 126]. Null geodesies are invariant w.r.t. 

conformal transformations, and Niederer's transformation, 

tanwt ^ X ujr^ 

I = , A = , b = s tan cur (6.1) 

u cos ut 2 

maps indeed every half oscillator period conformally onto the space-time which describes a 
free particle, 

dX'^ + 2dTdS = — I — (dx^ + 2dtds - u^r'^dt^) . (6.2) 
cos^ut 

This trick can not work for an anisotropic oscillator, though, otherwise the latter would also 
carry a full NH symmetry including rotation. 

An anisotropic oscillator is described by the metric [35] 

dx^ + 2dtds - {ulxj + ujlxl)dt'^. (6.3) 



Applying Niederer's transformation (6.1) i.e. 



arctanwT X ^1 u'^X^T , 
^= , X = , =, s = 5+ — (6.4) 

with some u then yields 

( dX^ + 2dTdS - -^d^^X'.dT' - -^d^^XldT' 

Now choosing either u = Ui or u = u)2 eliminates one, but not both oscillator terms, leaving 
us with 



df = 1— ( dX'' + 2dTdS — — ^^X'fdT' 



( dX^ + 2dTdS - -^^^LJ^xMtA . (6.5) 



= dX' + 2dTdS- . 

[where we should have put indices 1 or 2 on X, depending on our choice of u]. For both 
choices we get, hence, a maximally anisotropic "one-sided" "Hill-type" system, with Newton- 
Hooke symmetry — except in the isotropic case 

uJi = 002, (6.6) 

when both oscillator terms drop out, leaving us with a free system carrying its full Galilei 
symmetry. The latter can then be "re-imported" through the inverse of the Niederer trans- 



formation (6.1) to yield full Newton-Hooke symmetry. 

In conclusion, the "prototype system" is of the "Hill type" , with its rotation-less Newton- 
Hook symmetry — which, in the isotropic case, degenerates to a free particle with restored 
rotational symmetry. 
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7. CONCLUSION 



Souriau [S] attributes the center-of-mass decomposition of a free non-relativistic system 
to Galilei symmetry, more precisely, to an invariant Abelian subgroup of it, whose existence 
is rooted in turn in the cohomology of the Gahlei group [8j. Remarkably, it is this same 
cohomology which rules central extensions |9]. 

In this paper we performed an analogous study in the Landau problem, based on the 
Newton-Hooke group. The clue is that Newton-Hooke and Galilean symmetries are indeed 
"hiddenly the same" j3|, and have therefore identical cohomological structures [T^ . 

The intuitive content of Kohn's theorem, i.e., the relation between [Newton-Hooke] sym- 
metry and center-of-mass, is now clear : each particle, taken individually, would carry such 
a symmetry; Kohn's condition is precisely what is needed to extend this symmetry to the 
center-of-mass, which will hence represent the motion of all particles collectively. 

A method for finding approximate solutions of the 3-body problem of Celestial Mechanics, 
also used in Galactic Dynamics [27], is referred to as the Hill Problem. It was found recently 
[0] that the latter also has a symmetry reminiscent of the Newton-Hooke one, except for 
rotations, which are broken. 

At the technical level, the Hill Problem is a particular case of an anisotropic harmonic 
oscillator in an effective magnetic field. 

In this paper, we performed a similar study for a general anisotropic harmonic oscillator. 

All our investigations here have been purely classical. The decomposition of Newton- 
Hooke symmetry into Heisenberg algebras is, however, particularly useful for the quantum 
description, see [iHl [20] for details. 

Appendix A 

We find here the canonical form of the 4x4 antisymmetric nonsingular matrix u under- 
going the transformation 

u DujD'^, where D obeys D'^XD = X, (A.l) 



X being the matrix defined in eqn. (5.3). 



As it has been noted in the main text, X can be put into canonical form, which depends 
on the assumption concerning the eigenvalues of X. 
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Consider first X positive definite; then we can put X = I. As a result D is orthogonal and 



we have to find the canonical form of u under o(4) transformations (A.l). This resembles 



the problem of classifying the electromagnetic field configurations under the Lorentz group 
0(3, 1). Guided by this analogy, we define 



fi — <^Qi 



9i — \^jk'-^jk- 



(A.2) 



Note that fi and gt transform like vectors under SO (3) transformations acting on the last 
three coordinates. Moreover, det a; ~ (/ ■ gY, so that f ■ g 0, i.e., f 7^ 0, g and / is 
not perpendicular to g. 

Let us consider the rotation in the plane spanned by the 0-axis, and the axis which is 
orthogonal to it and defined by the unit vector n. The transformation rules under such a 
rotation read 



f\\ = f\\, /i = /_lcosv5 + (n X ^j_)sin(^. 



^11 = ^11' 



g± = g± cos If - {fix fx) sin (fi, 



(A.3) 



where || (_L) denotes the component parallel (orthogonal) to n. li f 1^ g we put 



n 



f ^9 
\fx9\ 



and 



sin2y9 



2\fxg\ 
P + 9 



(A.4) 



to achieve / || g. Then by S0{?>) rotation one gets further = VLi5i2, > 0, gi = —Q2Si2- 
Renumbering, if necessary, 1 o 3 (which is an 0(3) transformation) we let ^2 > 0. Due to 



definition (A.2), 



/ /i /2 /s \ 



/ \ 



in the form X 



-/i gs -g2 
-f2 -93 gi 
V-/3 92 -9i / 



. Put 



(A.5) 



1^2 
-fil 
\ -Q2 0/ 

Consider next the case of semidefinite X with one null eigenvalue. Then X can be put 




' d A 
B U 



D 



(A.6) 
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Eqns. (A.l) implies B = 0, U G 0(3), so D acquires the form D = \ \ , d 0. Then, 

U 



with Uij = eijkQk, 



DuD^ 







(A.7) 



Here Ug is an antisymmetric matrix, so it belongs to the algebra so(3). One can choose 
therefore U G SO (3) such that 



^ 1^2 ^ 



UUgU^ 



\ 





-^]2 



^2 > 0. 



Consider now the elements dfU'^ + AugU^ = dfU^ + AU'^UugU'^. Call 



(A.8) 



dfU^ = (/i,/2,/3), AU^ = (ai,a2,a3). 



(A.9) 



Then 



dfU"" + AU^UcOgU'^ = (/i,/2,/3) + (ai,a2,S3) 



^ 0^2^ 



-^2 



(/l, /2, /s) + (-^^203, 0, n2di). 



(A.IO) 



Knowing /, U and d one determines /i,2,3 and chooses ai^3 in such a way that 



rf/f/^ + AU'^UugU'' = (0, /2, 0), /2 ^ 0. 



(A.ll) 



By an appropriate choice of d we get < /2 = ^i, so (A.7) acquires the form (A. 5). 



Finally, let X have two positive, one negative and one zero eigenvalue. Without loosing 
generality, we put 



X 




G 



G = diag{—l, 1, 1) 



(A.12) 



With D of the form (A. 6) eqn. (A.l ) yields B = 0, f/ G 0(2, 1) ; DluD^ has the same form 
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Consider now UoOgU'^ . Again proceeding along the same lines as in the classification 
of electromagnetic field configurations, we find that UugU'^ can acquire three "canonical" 



forms: 



UUgU^ 



UUgU^ 



^ ^2 ^ 



-^]2 

^ 
A 
-A 







n2>o 



(A.13) 



A > O 



(A.14) 



UUgU^ 



\ 



(A.15) 



s 

-S -S 

If (A.13) holds the same reasoning as previously leads to eqn. (A. 5). In the second case 

/ n \ 



UJ 



-n 

A 

\ -A y 



n>o, A>o, 



(A.16) 



Finally, if A.15 holds 



/ n \ 

s 

s 
\ -s -s / 



(A.17) 



Appendix B 



We consider here the orbit method for the Lie algebra (5.1). The general element of the 
dual space can be written as 

hH + mM + ZiC- (B.l) 
Consider the coadjoint action of g = exp {iy^C,k)- It reads 



m = m. 



z[ = Zi + Ukiv'^m, 



(B.2) 



1 



h' = h + y'^uJkiXijZj + -y^y^uji^ujkjX^jm. 
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Assuming m 7^ and using the fact that u is invertible, we conclude that each orbit 
contains the points corresponding to = 0. The set of these points forms the coadjoint 
orbit of the stabihty subgroup of the relations Zi = 0. However, the latter is generated by 
M and H, so the coadjoint orbits are trivial. We conclude that Zi = define exactly one 
point on coadjoint orbit. Therefore, generating the whole orbit by the action of our group 
on that point we conclude that the orbit can be parametrized by the variables Zi and 

h = e + ^XijZiZj, (B.3) 

where e is the value of h at the point Zi = (internal energy). The basic Poisson bracket 
reads 

{zi,Zj] = UijiTi, (B.4) 

which completes the description. 

The additional symmetry generators can be dealt with in a similar way. 

In the case of two degrees of freedom and (semi) definite H it is convenient to identify 



the "physical" generators as described by eqns. (5.8) (i.e. to single out the boosts and 
momenta). In this basis the counterparts of dual coordinates Zi are denoted by pi and bi (cf. 



eqns. (5.11) and (5.12h). 
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